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(~| , Abstract 



We calculate the 0(as) gluon Bremsstrahlung corrections to the inclusive decay h — > 
57, involving the full operator basis Oi - Os- Confirming and extending earlier calculations 
of Ali and Greub, we give formulas for the total decay width as well as the perturbative 
photon spectrum, regarding the former as a necessary part of the forthcoming complete 
NLO analysis. We explore in detail the renormalization scale dependence of our results 
and find it considerably increased. 



1 Introduction 



If we define rare weak decays of hadrons as weak decays which are rare because they are loop- 
induced (as opposed to CKM-suppressed) , the analysis of these decays may shed hght on at 
least three important topics: (i) Since they can proceed even at leading order only through 
diagrams with a loop of virtual particles, these decays do test the standard model of weak 
interactions as quantum field theory, (ii) Since they are decays of hadrons, they are of course 
strongly affected by QCD effects. If, however, one naively calculated these corrections, due to 
the massiveness of the vector bosons in the loop one would have to cope with large logarithms 
In^^^, where Myy is the electroweak scale and nig the usual hadronic scale. These large 
logarithms spoil perturbation theory: one has to resum them invoking the powerful techniques 
of the renormalization group. Especially in the case b s'j, where the strong interaction is 
known to enhance the decay rate by a factor of 2-3, we have hence an ideal testing ground 
for those frequently employed concepts of partial resummation of the perturbative series, (iii) 
Since they are rare even in the standard model, they are sensitive to the effects of new physics: 
plainly spoken, new heavy particles running in the loop will contribute to the decay if their 
masses are not much larger than the electroweak scale, as expected for supersymmetry, left- 
right-symmetric models, technicolor and so on. Of course also the couplings of these particles 
must not be too small: so if we can find no positive sign of new physics, i.e. a deviation from 
the Standard Model (SM) prediction, some more or less sharp restrictions on the parameter 
space are all we can hope for. 

Among all rare weak decays b ^ s'-f takes a special role: it is of order G'j^^a (not Gp or 
G^a^), and accordingly the corresponding branching ratio is much larger than that of most 
other rare decays. In fact, within the field of B physics, 6 — S7 is the only one which has 
been measured experimentally, and this only recently. Following the first observation of the 
exclusive mode B —>■ K*^ in 1993 the CLEO collaboration reported by now measurements 
of the inclusive branching ratio, BR[i? Xg'-f] = (2.32±0.67) x 10~^ as well as of the photon 
spectrum. Here Xg denotes an arbitrary state of total strangeness —1, and experimentally some 
lower cutoff in the photon energy has to be imposed in order to exclude the tree-level channel 
b ccs'-f. As it is well known, the inclusive rate is of much more theoretical interest than the 
exclusive modes, because within the framework of Heavy Quark Effective Theory (HQET) one 
can show that it is given by the free quark decay model plus perturbative QCD, up to calculable 
corrections of order 1 / ml . The same is surprisingly true for suitable defined moments of the 
photon spectrum [Q. Recall that a continuous photon spectrum in 6 — > 57 arises by reason of 
two effects: perturbatively, by emission of gluon bremsstrahlung (the topic of our paper), which 
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results in a long tail of the spectrum below the kinematical endpoint, and non-perturbatively 
by the Fermi motion of the b quark inside the B meson, which leads to a symmetric smearing 
around the endpoint. Combining these two effects is a quite non-trivial task, primarily because 
one has to make a consistent separation between the perturbative and the non-perturbative 
region. The analysis of the spectrum has been the subject of some recent papers ||^, ^, |^ ^, 
and it will thus not be the central point of our work; instead we will give complete results for 
the bremsstrahlung contribution to the total decay width. To our knowledge, there exists only 
one calculation of Ali and Greub on this subject, which for a long time was in need of an 
independent check. In view of the increasing desire for precise predictions of BR[i? —* Xs'~f] we 
have performed this check, confirming and extending Ref. f^. 

Our paper is organized as follows: In section 2 we give the effective Hamiltonian which we 
used in our calculations, and discuss very briefly the structure of a complete Next-to-Leading- 
Order (NLO) analysis and the place of this paper within that greater task. In section 3 we give 
the amplitude b s'jg using this Hamiltonian, and we also calculate those 0{as) corrections 
to the amplitude b ^ s'j which are needed for cancellation of infrared divergences. In section 4 
we compute photon spectrum and total decay width from this amplitude, giving some details 
about the phase space integration in Z) = 4 — 2e dimensions. Subsequently, in section 5 we 
analyze these results numerically, exploring their range of applicability and discussing their 
dependence on the renormalization scale as well as on the other input parameters. The paper 
closes with a short summary and outlook in section 6. 



2 The effective Hamiltonian 

In order to make use of the renormalization group techniques for calculation of short- distance 
QCD effects, we work within the framework of an effective five quark theory where the W boson 
and the top quark have been removed as explicit dynamical degrees of freedom. Neglecting 
contributions with smaller CKM parameters (iKife^Tsl/lV^feV^sl < 0.02), the relevant Hamiltonian 
for the process b — > 57, sjg is in leading order of the operator product expansion given by 

Heff = -^V,,V*J:C^i^^n. (1) 

Here, Gi;' is the Fermi coupling constant, Vij are the CKM matrix elements, Ci{n) denote the 
Wilson coefficients evaluated at the scale /x, and the operator basis reads 

O2 = {cLal^hLa){sLI3lfiCLf}), 
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q=u,d,s,c,b 
q=u,d,s,c,b 

q=u,d,s,c,b 
q=u,d,s,c,b 

6 

Os = -^fspa^^{m,R + m,L)T|„6,G;:„ (2) 

with e and Qs denoting the electromagnetic and strong couphng, respectively, and the pro- 
jection operators R,L = |(1 ± 75). Some caution is necessary if we apply this result, which 
was originally obtained only for the process b —>■ 57, to the process b s'yg. This is, because 
potentially new diagrams like the one shown in Fig. 1 could lead to new operators containing 

b ^ 



Figure 1: A standard model Feynman diagram contributing to the decay b — > s^g. Diagrams of this 
type lead to low energy effective operators which are not contained in the Hamiltonian (1), which was 
originally derived for the process 6 — > 57 only. However, these operators are of dimension 8 or higher 
and hence are suppressed by powers of ml/m^r. 

two quark fields, the photon and the gluon field simultaneously. However, a detailed analysis 
following the lines of |jlO[, i-e. using gauge invariance and the equations of motion, shows that 
all these operators must necessarily be of dimension 8 and hence are suppressed by powers of 

The inclusive decay rate was calculated in the approximation of leading logarithms, i.e. 
holding terms of order (a^ In M]y / mi,)"' , n = 0, 00. An exhaustive discussion of these results 
can be found in Ref. In this reference, there is also analyzed the general structure of next- 
to-leading calculations, which hold terms of order as{asln.Mw/mh)"',n = 0, ...,00. The main 
points we would like to stress here are: (i) The dominant theoretical uncertainty in the leading 
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order result arises from the renormalization scale dependence of the Wilson coefficient C-j{^). 



It amounts to about ±25% [|T2| , and it can hopefully be reduced by the full NLO calculation to 
about ±3%. (ii) The full NLO calculation requires three steps: 0(a<() matching of the effective 
theory to the SM matrix elements at /i = Mw which gives Cj(/i = M^y), evolution of these 
Wilson coefficients down to yU ~ mf, with the help of the NLO anomalous dimension matrix, 
calculation of 0{as) corrections to the matrix elements of the effective theory at the low scale. 



From these three steps the first has been tackled in Ref. |T3[ . Certain parts of the second have 
also been calculated, see Ref. for a summary of the present state of the art. Unfortunately, 
the most difficult (and presumably fairly significant) ones, involving finite parts of two-loop 
and divergent parts of three-loop diagrams, are not yet known. The last step splits into two 
pieces: 0{as) corrections to the b ^ s'j amplitude, and the bremsstrahlung amplitude b S'jg 
to 0{as). The latter has to be considered, because the final state s+gluon contributes to the 
measured state Xg as well as a single strange quark. In the following sections, we calculate 
these bremsstrahlung corrections using the full operator basis (2). We would like to emphasize 
that the contributions of the penguin operators O3 - Oq, which were neglected in 0, should 
not be considered numerically negligible to begin with although their Wilson coefficients are 
small: in leading order, their effect on the decay width amounts to about 15% (in the NDR 
schemeQ), hence a consistent calculation certainly should include all their contributions beyond 
leading order, too. 

We now proceed to give a somewhat detailed calculation of the above mentioned brems- 
strahlung contributions, closing with some numerical results at the end of section 5. 



3 The amplitude b s^g 

Using the effective Hamiltonian (1) we will calculate in this section the complete transition 
amplitude for the decay 

where a, j3 and a denote colour indices and e and 77 are the four dimensional polarization vectors 
of photon and gluon, respectively. The amplitude will have the form 

M"--^ = E M, = T|„\/^ a(/i)n(p')T.n(p), (3) 

i=l i=l 

Unless noted otherwise, by scheme we mean throughout this paper the regularization scheme used for 
the treatment of 75 in D dimensions. NDR stands for Naive Dimensional Regularization, and HV for the 
't Hooft-Veltman scheme. 
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Figure 2: The tree- level Feynman diagrams which contribute within the effective theory (1) to the 
decay b ^ s^g. 

where we have extracted the couplings in a common factor 

denotes as usually the Gell-Mann matrices, Tj are some Dirac structures depending on p,q,r, 
e and rj, and Ci are the so-called effective Wilson coefficients [|ll| which are linear combinations 
of the usual Wilson coefficients (|l]), defined (in the NDR scheme) by 

C'i...6 = Ci...6; Cj = C7 — -C5 — Cq, Cs = Cs + C5. (5) 

In the HV scheme, Ci = Ci. The reason for introducing these quantities will be explained after 
calculating the matrix elements of the penguin operators O3 — Oq in Sec. 3.3. 

Throughout this section, we work (where necessary) in D = 4 — 2e dimensions in order to 
regularize UV divergences, and we treat 75 in D dimensions according to the NDR scheme. 

3.1 Magnetic penguins: O7 and Og 

The magnetic penguin operators O7 and Og contribute only through the tree-level diagrams of 
Fig. 2 which yield trivially 



^ 2 

rri Qd 



-{mbR + msL)(^li{i) -f + mb)^ + ^^Hi - ^ + ms)^^{mbR + m^L) 



-2pr 2p'r ' 

-{rribR + msL)fi]{^ - i + rrib)^ + -^i^{1> -f + ms)f^{mbR + nisL) 



(6) 
, (7) 



-2pq 2p'q 
with Qd = —1/3. 

Now these parts of the bremsstrahlung amplitude stemming from the magnetic penguins 
are affected by infrared divergences. That is, after phase space integration IM7P and iMgp 
yield no finite decay rate. Therefore the first step is to regularize these divergences: for the 
sake of simplicity, we choose dimensional regularization [|T^ so that the phase space integration 
is to be done in Z) = 4 — 2e dimensions and therefore becomes a little bit involved. The second 
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step is to eliminate the divergences. According to the well-known Bloch-Nordsieck mechanism 



16| , this can be achieved by considering the appropriate diagrams with virtual instead of real 



gluons. Since the infrared divergences of Og occur only at the low energy endpoint of the photon 
spectrum which is experimentally not accessible, they are only of minor interest and need not 
to be considered here. We are left with the one diagram of Fig. 3f which yields (together with 
the 0(a°) contribution) in the Feynman gauge and before renormalization: 

= ^^V,,V:fiM5^,u{p')Ti%,Mp). (8) 



with 



^r = ^(4-)T(l + e)l±^ 



2 

mt , , 

P = —2- 9 

ml 

Note that this amplitude is UV finite: the 1/e poles correspond to infrared singularities. If 
we wish to calculate the transition amplitude, we also must make allowance for diagrams with 
quark self energy parts. To this end, remember that according to the LSZ reduction formula the 
transition amplitude is given by the corresponding truncated Green function times the residua 
of the propagators of the external particles to the same order in a^. Now let us employ the 
on-shell renormalization condition 

^l.™-0 (10) 
for determining the finite parts of the wave function renormalization constant Z^, which yields 



Z4m) = l-^(--37 + 31ni;^ + 4 (11) 



in Feynman gauge (here the singularity contains an 1/e UV and an 2/e IR contribution). With 
this choice, the residuum of the renormalized quark propagator equals 1 (to one loop), and 
the transition amplitude is simply given by the renormalized truncated Green function, i.e. by 
the amplitude (^ given above after including renormalization. To achieve the latter, we have 
to add the counterterm amplitude {Z-j'jZmZ]L'^{mi,)Z]l'^{ms) — l)Mr^, where we choose the 



operator and mass renormalization constants Zjj and Z^ according to the MS scheme |17 



At this place, employing the MS scheme is of course necessary, because the whole calculation 
of the Wilson coefficients has been done in this scheme. The result of this renormalization 
procedure amounts to replacing T^^^irt in (H) by 

Ti%t = irnt{fi)R + ms{fi)m{l + irf )), (12) 
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(a) 



(b) 



(c) 






(d) 



(e) 



(f) 



Figure 3: Examples of one-loop Fcynman diagrams which contribute within the effective theory (1) 
to the decay b s'jg. (a), (b) and (c) depict insertions of the current-current operators, whereas (d) 
and (e) illustrate the two distinct types of insertions of penguin operators, (f) displays the virtual 
gluon correction to the b ^ S'j matrix element of the magnetic penguin operator O7. This graph is 
needed for cancellation of infrared divergences at the high energy endpoint of the photon spectrum in 
b sjg. 



with 



Kf>= |(4.)T(l + e)|-i 



P 



mz 



l±p 
1-p 



lnp + 2 



1 



mt 



2 2 ^ 



m 



P^ 



mt 



+ -ln-f + -ln-^-4 



P^ 



(13) 



and mq(iji) the running quark masses in the MS scheme. Again, this is UV finite. 



3.2 Current-current operators: Oi and O2 



The current-current operators contribute through diagrams of the type shown in Fig. 3a. Note 
that diagrams as Fig. 3b in which only one particle is radiated from the loop vanish on-shell, as 
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one may easily verify through exphcit calculation. From chirality, the same is true for diagrams 
as Fig. 3c in which there is no particle at all radiated from the loop. 

Since their colour structure entails a factor tr T°, all diagrams with an insertion of Oi vanish: 

Ti = 0. (14) 

Using some easily calculable one-loop three-point functions, one obtains for the diagrams in- 
volving O2 



m 



Q^K W2, (15) 



with 



ifj = 1^ [{rnm - (t'-Wf - (crXwXrf - ;')] + Mi - f)- ('-'iM -f) + 2(w)4 i, (w) 
= i + (17) 

(18) 

Qu = 2/3 and the function Qq{s) defined by 

Qq{s)= —ln{l-sx + sx^-i6), (19) 



/O X 

with 5 positive infinitesimal. The integration in ( |T9|) can be performed analytically, yielding 
Qo{s) as given in eq. (9) of Ref. 0. Note that all UV divergences cancel to yield an UV-finit 
result. Note also that (H), (|^) and ([T5|) are in complete agreement with Ref. after some 



rearrangements in the Dirac structure using the Dirac equation. 
3.3 Penguin operators: O3 — Oq 

Evaluating diagrams containing penguin operators, one must take into account two types of 
possible insertions: the somewhat natural one of Fig. 3d, but also the one that leads to diagrams 
like Fig. 3e. The latter stems from terms in the second current of the operator containing the b 
or the s quark. But as it turns out, only one insertion for every operator survives due to colour 
structure. 

In a calculation analogous to the previous one and without making a Fierz transformation 
we obtain: 

T3 = Q4tib + ^W2, (20) 

T4 = {l + Qd[^b + ^s]+Qu'ic]w2, (21) 



6 

T5 = Qd{mhW^R + msW^L), (22) 
Te = -T4, (23) 
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where we have used the abbreviation k, 




. W2, ti{s), and s are defined in ([T6|) , ([17|) 



, respectively, and 

(24) 

There is one more point to be considered: in the case of the penguin operators O5 and 
Og, there are also non- vanishing diagrams in which at least one photon (or gluon) is radiated 
from an external leg. But as in 6 — 57, sg the resulting amplitudes are proportional to either 
T7 or Tg. Thus the same redefinition of the Wilson coefficients as in 6 — s'y,sg absorbs 
these contributions into the tree- level amplitudes T7, Tg. This is why we have introduced the 
quantities Ci as defined in (|^)0. 



+ 4 



(r/e) - 



qr 



4 Spectrum and decay rate 

We now proceed to calculate some physical observables from the amplitude M'"'^'"^ combined 
with 



4.1 Squared and summed amplitude 

Summing the squared amplitude |]Vf^^'^'^*p over spins, polarizations and colours is a straight- 
forward procedure. First of all, we can write 

|^6rems|2^^1 J2 \M'^ems^' = ^mt\V\' Qif^)C,ifi)M,^, (25) 

spin, poL, col. i.j—l...s 

j>i 

where the factor 1/6 stems from averaging over spin and colour of the incoming b quark and 
V is given in (^. Introducing the kinematical variables 

2qr 2pr 2pq 
s = — 77, t = — 5-, u = — 5- (26) 

and again using the abbreviation s ) , the contributions from the various operators 

can be expressed in the limit = as 

Strictly speaking, these diagrams contribute only in the NDR scheme; in the HV scheme they vanish and 
therefore in the latter Ci — Ci. This difference is compensated by a corresponding scheme dependence of the 
leading order anomalous dimension matrix. 
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M22 
M23 
M24 
M25 
M27 
M33 
M34 
M35 
M37 
M44 



M47 

M57 
M78 



9 

-^Re 
9 



8, 



-M26 = -Re [kI{2kc - Kb)] (1 - s), 
9 

-^Re[<(l-2«6)](l-s)s, 



28 



-Re{Kc)s, 



1,1 ,9/ 



9'2 
-M36 

-Re 
9 



--Re 
9 



{2 + '^Di'^'^c- i^b) 



(1 



-3M38 = ^Re(^ + Kb)s, 



66 = ^|2/tc - /tfcPll - S), 



M45 = -M56 = --Re [{2kc - Kb){l - 2kI)] (1 - s)s, 
9 



-M67 = -3M48 = 3M68 = --Re(2fi:c - Kb)s, 
lel^bl^ + 1(1 - '2Kb)s + AKbl'^] (1 - s), 

3M58 



--Re(fi;6)s, 



2 \ 



(27) 



For rus ^ the corresponding expressions become quite lengthy. The complete formulas for 
this case can be found in the appendix. 

There are still two contributions missing. M77 is the one which will by far dominate the 
photon spectrum (at least near the endpoint) , so in this case it may be interesting to know the 
full formula (i.e. 7^ 0), which reads 



M77 = (1 + p) {(1 - p)m(^) - 2Mif} 



(1) 



77 



2u(u-2)- 2u-l^ 
1 + M H 7 — + -P 



1-u 



^77 - p 



l-^t + 



u 



{1-u 



r,\2 



1-U 



(28) 
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where we have introduced the rescaled kinematic variables i = t/{l — p), u = u/{l — p), and 
p = m?g/m^ from (^). In writing (pSf), we have spht M77 into an infrared safe part M^y'' and 
an infrared divergent part M^'^ . From the amphtude it should then be obvious that the 
quantity Mss/Qd is given by the above expression for M77 if one interchanges u with i. 



4.2 The infrared finite parts 

As stated above, in the region of experimental interest all contributions Mij to the squared 
amplitude yield a IR finite spectrum and decay rate, with exception of M77. For those finite 
contributions, we therefore can do all phase space integrations in 4 dimensions. Note that from 
conservation of energy-momentum we have the relation u + t — s = 1 — p and the kinematic 
boundaries m G [0, 1 — p], t G [1 — u — p, 1 — j^]- Furthermore, in the rest frame of the 
decaying b quark the variables u and t become the energy of the emitted photon and gluon, 
respectively, measured in units of mb/2. Splitting the spectrum into an infrared finite and an 
infrared singular part, 

du du du ' 

one therefore obtains trivially for the former 

" E ^ = E A(p)A(/^) / rftM.,(t,«), (30) 



du u,.^ u/i 



with 

^ \VttV:fGla 



ml (31) 



327r4 

Regrettably, due to the presence of the functions n{s) in the expressions (^) for Mij the 
integrations in ( ^9l) cannot be done analytically, so we postpone the numerical analysis of (|29|) 
to sections 5.1 (spectrum) and 5.2 (decay width). 



4.3 The infrared divergent part 

Integration of IM7P yields a photon spectrum which contains a non-integrable singularity as the 
photon energy reaches the kinematical endpoint, —>■ ^(1 — p). To regularize this divergence, 
we will do the whole phase space integration in D = 4 — 2e dimensions. The singularity will 
then manifest itself in a part of the decay width F^*"^*"^ = J dV'^l^"'' proportional to 1/e. This 
singularity will cancel if we add the equivalent singular contribution F^*''* stemming from the 
virtual gluons, leaving us with an unambiguous finite result F7 = Py^*^"^* + F^*^*. Let us now 
calculate F7. 
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In D Dimensions, the Lorentz Invariant Phase Space (LIPS) is given by 



(32) 



(27r)^-i2po' 

and accordingly the differential decay width for an n-body decay with transition amplitude M 
reads 

1 f n n 

(33) 



= — (27r)^5(^) \p-Y.P^j \M\'X{d^{p.). 



1=1 

Specializing this to the case h s'yg and performing the delta function as well as the angular 
integrations (note that our squared and summed amplitude |M^''^™"^|| ( [25| ) only depends on the 
energies of the particles in the final state), one obtains 

2e 



dr 



brems 
7 



2a. 



-4e iin^J.^ 



37r"" r(2-2e) 
where we have defined the function 



-2e 



tu{l - z{t, uy) MTjit, u) dt du, (34) 



u + t-1 
ui 



(35) 



which is actually the cosine of the angle between q and f in the rest frame of the decaying h 
quark. The crucial point to be observed now is that even in the first integration over i (which 
yields no divergences) one has to hold terms of order e, because they may multiply later on 
with 1/e poles of the u - integration and so give finite contributions to the total decay width. 

After carrying out the integration over t, we obtain the following expression for the photon 
spectrum to 0(e): 



^^(i + p)(i-p)3c,^(.)|iroa 



5(1) 



2^(2) 



l-U 



(36) 



with 



2e 



C(2) 



-U 



^2t 



g-p)-" (^) 

r(2 - 2e) 

1,^ .u{l-u){2u-l) ,m(2m2-5m-1) ,^ ,^ 

2 (1 — uy 2 1 ^- 

5(2) + eS^^ 



u 



b ' 



l-U 



1 + p 

u H ln(l — u) 

1- p 



S, 



(2) 



^ [2 + ln(l-M)]M-2^^ ^ln(l-M)+ ^"^^ 



u 



P 



P 



-ln2(l -m) -2 Li2(M) 



(37) 
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p = ml /ml and Li2(x) = — /q y ln(l — xt) is the Spence function (or dilogarithm) . From this, 
we can now accomphsh our final goal, the integration over u. Since our results agree entirely 
with (29) and (30) of Ref. 0, we will not give them here explicitly. 

Then one has to go through the same steps with M^*^*: square it, sum over spins and 
polarizations, do the phase space integration in D dimensions. The result will contain a singular 
term that cancels the one of r^^^'"*. In fact these calculations are much simpler than the one 
sketched above, because we only have to struggle with a two-body decay. After some work, we 
arrive at the following final resultQ for the total decay width via O7: 



{31 



with 



and 



mb[mb) 



(39) 



ST7 




4 2 - 7p + 2p2 p(7 - 8p + 5p2 



(1 - Pf 



Inp - 41n(l - p) 



4 Li2(p) - -TT^ + 2 In p ln(l - p) 



+ 41n 



mt 



p^ 



r 



(0) 



(40) 



In deriving (pq) - (^Of) we have expanded the running quark mass mb{p), which occurs from 

17 , diuuina /X — mb, 



T2I) in the matrix element M"***, around p 



mb{p) = mb{mb) 1 + 



2as{mb) 



TC 



In^ 



(41) 



and we have also made use of the equivalence of mb{mb) and the squared pole mass ml up to 
corrections of order as{mb)- In the limit = (|38D looks much nicer: 

2a, /8 - 271^ 



c,{py 



37r 



+ 41n 



mt 



p- 



mb[mb) 
, rUb 



0- 



(42) 



We now hasten to undertake a critical discussion of all these results. 



■^^The attentive reader will perhaps feel uneasy about the fact that in order to obtain the following result, we 
apparently have calculated F"'"* with the effective Wilson coefficient Cy instead of C7. Of course the full two- 
loop matrix elements of O5 and Oe are not proportional to M™"'^ . Nevertheless, due to cancellation of infrared 
divergences they must contain at least a part proportional to Mj"'*. The other parts have to be considered in 
the full calculation of two-loop b 57 matrix elements, which is not the aim of this paper. 
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5 Numerical analysis 



5.1 Photon spectrum 

If we neglect all non-perturbative effects, what does the photon spectrum due to bremsstrahlung 
actually look like? The answer depends on the experimental setup: Let us assume the resolution 
of the photon detector is so that in the endpoint region we cannot discriminate photons 

2 2 

within the energy interval E™"^ — < < E^"^, with E^"^ = ""^"^^ the maximum 

energy of the emitted photon. If we now define the function 

Ttotiuo) = Tr - / du + / -j^du, (43) 

Jo du Juo du 

which is the total decay width including all photons with an energy > ^Uq, then the photon 
energy spectrum is given by 



dr'' 



dTA 

du 



du ^ n < 1 - p - A 

ir^ot (1 - p - A) , l-p-A<u<l-p (44) 
0, u> 1- p 



Here, dV^^'^^^/du is the perturbative bremsstrahlung spectrum which we calculated in Sec. 
4, cf. eqs. (^9]), ( ^OD and (|36D . We have plotted the resulting photon spectrum in Fig. 4, where 
we have used = 200 MeV, A = 25 MeV and the Wilson coefficients Ci{mf,) as given and 
discussed in the following section. Though we will not dwell upon any further details, a few 
comments may be useful: 

First of all, if we consider the limit A ^ 0, up to what minimal resolution Acrit the spectrum 
(|4^) can be trusted? Clearly, this is a question of the reliability of perturbation theory. We can 
estimate Acrit if we consider the case where the pure 0{as) bremsstrahlung contributions from 
the region u < 1 — p — A become of the same order of magnitude as the contributions from the 
peak in the region 1 — p — A<m<1 — p: 



1 p ^crit ^Ybrems 



I 



du 



du ~ Ttotil- P- Acrit) ■ (45) 



This is, because if the l.h.s. of (^) becomes greater than the r.h.s., then the 0{as) brems- 
strahlung corrections exceed the two-body decay peak so that perturbation theory is no longer 
valid. If one wishes to explore the case A < Acrit, one will have to resort to some Sudakov-type 
partial resummation of the most singular terms ( "exponentiation" ) in order to obtain a reliable 
spectrum. However, numerically we get from (^5]), with rris = 200 MeV, Acrit — 0.015, so that 
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Figure 4: Two different views on the B Xg'j plioton energy spectrum due to gluon bremsstrahlung, 
eq. (|4^). The spectrum is given in units of |C'7(mb)pro, the leading order decay width with rus = 0. 
In both plots, the solid curve corresponds to the full operator basis and nis = 200 MeV. (a) shows the 
region 0.6 < u < 0.95, where u is the energy of the emitted photon in units of mi,/ 2. The dotted curve 
represents the contribution from the operator O7 alone, which is obviously the dominant one. (b) 
takes a closer look at the endpoint region, where the spectrum is sensitive to the mass of the strange 
quark. Here, the dashed curve depicts the case = 500 MeV. 

an infrared improvement remains unquestionably necessary only in the last, say, 30-40 MeV 
near the kinematical endpoint. This is well below any experimental resolution, and moreover 
it will presumably be of no importance for the lower moments of the spectrum ||^. Therefore, 
we suggest to take ( ^Ij ) with A ~ 0.015 as the perturbative photon spectrum and as the input 
for the inclusion of non-perturbative effects. 

From Fig. 4a one sees that the spectrum is completely dominated by the contribution 
stemming from the Operator 0^ (solid vs. dotted line). We can express this observation in a 
more quantitative way, if we define the (normalized) moments of ( ^^ as 

M^iuo, A) = — — - / m"-^ du, n G N. (46) 

LtotyUojJuo du 

For instance, we obtain thereby for the first moment Mi (0.6, 0.015) = 0.961 if we include 
contributions from all operators, and Mi(0.6, 0.015) = 0.958 with Cj(/i) = for i 7^ 7. These 
dominance of O7 is a quite welcome feature, because dTj^^^^'^/du is the only contribution which 
we can calculate analytically (without any numerical integration). 

The moments (|46|) are also useful to illustrate the dependence of the spectrum on the mass 
of the strange quark. Varying between 200 and 500 MeV (solid vs. dashed line in Fig. 4b), 
Mi(0.6, 0.015) changes from 0.961 to 0.952. So the moments are decreasedhj increasing values 
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of rUs, in contrary to the result found in Ref. R. This is primarily due to our slightly different 



normalization of the moments, eq. (0), which we have chosen because the so-defined moments 
can in principle be extracted from experiment without knowledge of rris- On the other hand, 
it is also due to our different treatment of the infrared divergent endpoint region: whereas 
Kapustin and Ligeti trusted the perturbative bremsstrahlung spectrum up to arbitrary small 
A, with (44) we rely on first order perturbation theory only where it is clearly valid. Since the 



m^-dependence of the spectrum manifests itself mainly in the endpoint region, we regard our 
method as the more appropriate one. 

5.2 Decay width 

Taking into account the discussion of Sec. 2 it should be perfectly clear that from our cal- 
culations we cannot present any definite numbers for the total decay width (or, equivalently, 
the branching ratio). What we can do is, however, to list the complete bremsstrahlung con- 
tributions to the decay rate. When we will know all the remaining parts of the NLO analysis, 
especially the two-loop 6 — 57 matrix elements of the four-quark operators, we may then 
simply put together all these numbers to get some final answer. 

We will give our numerical results in units of \Ci{mi,)\'^rQ, which is the leading order decay 
width if one neglects the mass of the strange quark. Therefore, our only free parameter is 
p = ml /ml. Additionally, we have to assume some value for i?™", the minimum photon energy 
we can see with the detector. For some choices of p and -E™™, the resulting bremsstrahlung 
contributions to the total decay width calculated from (|30| ) and ( |38D are shown in Table 1. 

In our analysis we have used the leading order effective Wilson coefficients C*j(/i), see e.g. 



11 1 for explicit formulas, as well as the two- loop running coupling as(/i), and we set p = mi, 
with mh = 4.8 GeV being the pole mass of the b quark. The resulting numerical values for the 
Wilson coefficients Ci{mh) are displayed in Table 2, where they have been computed with the 
parameters Mw = 80.2 GeV, mt{mt) = 170 GeV and a,(Mz) = 0.117. 
Let us now discuss some interesting features of Table 1: 

(i) Contributions from different operators 

The dominant contribution to ST tot = Ttot ~ Ft-^'' comes from the operator O7 alone, all other 
contributions amount to only 15% of this one. However, this statement has to be taken with a 
grain of salt, because O7 is the only operator where we have included the virtual corrections. 
Contributions from the penguin operators O3 - Oq (last but one column) are completely negligi- 
ble as far as they are not included in the definitions of the effective Wilson coefficients. Whereas 
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■p 

J- 22 


■p 

J- 27 


■p 

J- 28 


p 

i 78 


p 

J- 88 


all others 


J- tot 




1.92 GeV 


66.92 


0.94 


-0.09 


0.03 


0.16 


0.02 


0.00 


67.98 


0.01 


1.44 GeV 


68.83 


1.85 


-0.27 


0.08 


0.32 


0.09 


0.01 


70.91 




GeV 


69.16 


2.48 


-0.48 


0.15 


0.48 




0.02 






1.92 GeV 


67.54 


1.14 


-0.08 


0.02 


0.20 


0.06 


0.01 


68.89 


0.002 


1.44 GeV 


69.09 


2.12 


-0.27 


0.06 


0.37 


0.18 


0.01 


71.52 




GeV 


69.40 


2.80 


-0.48 


0.10 


0.54 




0.02 





Table 1: Different contributions to the total decay width B Xg^ due to gluon bremsstrahlung, in 
percent of 1(77(7716)1^ Fq, the leading order decay width with rris = 0, as a function of the parameters 
p = ml /ml and i?^*". We have used the leading order effective Wilson coefficients from Table 2, 
mc = 1.4 GeV, mb = 4.8 GeV and rni,{mi)) = 4.4 GeV. The last column shows the sum of all the 
different contributions. 





C2 






c. 








-0.242 


1.104 


0.011 


-0.025 


0.007 


-0.030 


-0.309 


-0.147 



Table 2: Effective Wilson coefficients Ci{fi) at the scale /x = mb = 4.8 GeV, in the approximation 
of leading logarithms. These coefficients were evaluated with the two-loop as(/x) and the parameters 
Mw = 80.2 GeV, mtimt) = 170 GeV, as(Mz) = 0.117. 



the virtual corrections lower the decay rate significantly, bremsstrahlung tends to enhance it 
slightly. 

(a) Dependence on nig and E'™" 

Since <^ m^, from the outset we did not expect the mass of the strange quark to make any 
sizeable effect in Ttot-, and we used rris therefore mainly as regulator of coUinear divergences. 
Indeed it turns out that varying between and 500 MeV leads to a change of only 0.9% 
in Ttot- Here, we should note in particular that r7 = Fy''^™'^ -|- F"''* and all bremsstrahlung 
contributions Fj^ with ij ^ 77 are well defined in the limit — 0, except for Fgg. This 
divergent behaviour of Fgg has recently been discussed in Ref. [§], but since, from Table 1, F88 
is in the relevant part of the spectrum and for sensible values of of no numerical importance, 
this observation is of no great relevance here. 

The dependence on iv;™"- can be interpolated from Table 1. Note that values ii^™" < 
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1.4 GeV are not realistic. Above that, if one is interested in the hmit i^™" 0, one has to 
deal with the infrared (soft photon) divergence due to Os, cf. Ref. P, §]. 

(in) Dependence on the renormalization scale /i 

As mentioned in Sec. 2, an important offspring of the complete NLO calculation should be 
a considerably reduced renormalization scale dependence and thus a considerably improved 
theoretical uncertainty. From theory we know that in order to perform the correct resummation 
of large logarithms the renormalization scale ^ should be fixed to a value of order rrib, but that 
on the other hand it is arbitrary provided that In ^ remains a small parameter. Thus this 
uncertainty can be estimated if one varies ^ between, say, m;,/2 and 2m;,. Performing this 
exercise with the expression (^2]) for Fy, one obtains a startling large relative change of I'iggy". 
This is an even larger //-dependence than the ^^2^^ pure leading order result. What 

lesson can we learn from this peculiar effect? Of course a very elementary one: namely, that 
it is a mistake to expect any improved predictive power from partial NLO calculations. Let us 
explain this in detail. 

The //-dependence of the Wilson coefficients Ci{fi) is governed by the renormalization group 
equation [0 

/^^C'i(/i) = 7,i(«,)Q(/i), (47) 
where the NLO anomalous dimension matrix is given by 

7(a.)=¥°^^ + 7«f^V. (48) 

Here, 7*-°'' represents the leading order (LO) contribution, and 7*-^'' the NLO correction to it. 
7*^°) has been calculated by various authors and is given e.g. in Ref. |Tl|, and the calculation 



of 7(1) is not yet completed. From (0), ( ^HD we can infer that the dominant /t-dependence of 
6*7 (/t) near /t = is given by 

C,{^) = Crim,) + ^4^71? A(m,) In + O (a'. In A, «2 1^2 A) . (49) 

The second term in (|49|) is exactly the one that after a complete NLO calculation should be 
cancelled by the /i-dependence of the b ^ s-y matrix elements {Oi{fi)), since such a calculation 
should be free of renormalization scale uncertainties at 0(0;^). From these general grounds we 
therefore conclude that the transition amplitude 6 — > 37 at the NLO-level, 

M(h - .7) = -^V,,V,:5:a(/i)^^«(4(/x)), (50) 
v2 i=i 
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can be written as 



An 



-jPCiirrib) In — + ridinib) 



rUb 



{O7) 



tree 



(51) 



where rj are functions of the external momenta, spins and polarizations that contain no /i- 



dependence, and the last factor denotes the tree- level matrix element of O7. In (|5TD , the first 
summand C-j{^) has to be taken at the NLO level, whereas for the other Wilson coefficients 
LO is sufficient, since they are multiplied already with a^. In order to obtain a prediction for 
the /x-dependence of the quantity Vj which we calculated explicitly in this paper, cf. eq. (^2]), 
we square d^), set C'j(yu) = for i 7^ 7 and get 



1 + 



47r 



(0)i , ~ 

777 In ^ + r7 



(52) 



where f-j is an unpredicted number. With 777'' = y from |TT|, we indeed arrive at the same 
coefficient of In ^ as found in eq. (|42|) , and cancellation with the /i-dependence of Cj^ii) should 
work. Obviously, it didn't. This is for the following reason, which illustrates once more the 
importance of operator mixing under renormalization: since we neglected all the 6 — > 37 matrix 
elements (Oi(/i)) with « 7^ 7, only that /i-dependence of Cj{^) which is proportional to 777"*, cf. 
eq. 



, could be cancelled within our calculation. But 777'' C7{mb) < and J2 'la^Ci{mb) > 



\lri C7{mb)\, so actually C'j{ij) < is a monotonous growing function of /i. Consequently, 
by cancellation of the /x-dependence due to the anomalous dimension 777'', this monotonous 
growth is even more enlarged instead of being reduced. Such a reduction will only take place 
after including the 0{as) corrections to all b ^ s'j matrix elements {Oi{fi)). 



6 Summary 

In summary, we have calculated in this paper the complete 0{as) gluon bremsstrahlung correc- 
tions to the decay B X^'y. We confirmed all the results of Ali and Greub, and additionally we 
gave in Table 1 a compilation of the complete bremsstrahlung contributions to the total decay 
width. We also presented formulas for the photon energy spectrum and made a simple sugges- 
tion how to use this bremsstrahlung spectrum as an input for the inclusion of non-perturbative 
effects. We found that in both cases contributions from the QCD penguin operators O3 - Oq 
are entirely negligible, except for those parts of them which can be taken into account by the 
use of the effective Wilson coefficients (|^) instead of the usual ones. We investigated the consid- 
erably increased renormalization scale dependence of our results and showed that a reduction 
of this large theoretical error will not be possible unless all b ^ s'j matrix elements {Oi{^)) will 



19 



have been calculated up to 0{as). Very recently, the first calculation of such two-loop matrix 
elements was indeed performed [|l9l; as it should be, the results allow for a drastically reduced 



renormalization scale error. 

Finally, we would like to emphasize the urgent demand for the calculation of the still out- 
standing parts of the complete NLO analysis. Without these calculations, a consistently im- 
proved prediction for the branching ratio will not be possible, so we do not include such a 
prediction in our paper. We hope that these improvements in theory will be forthcoming soon, 
and that it will be possible to compare them with even more improved inclusive measurements, 
too. 
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Appendix 

Below we explain how to obtain from eqs. ( P7| ) the complete expressions (i.e. the case 7^ 0) 
for the various contributions Mij to the squared and summed bremsstrahlung amplitude (^Sf). 
First of all, in every factor 1/2 + Hf, and — i^b in ( pTj) replace by Kf, + ~ 1/2. 
Additionally, if both i and j G {2, 3, 4, 6}: make in the formulas of (^) the replacement 

{l-s)^{l-pf-{l + p)s. (53) 

If 2 7^ 5 and j = 7: make in the formulas of (p7|) the replacement 



If z 7^ 5, 7 and j = 8: make in the formulas of ( pTj) the replacement 

s-.il + p)s + -^^. (55) 
[S — u)u 

Finally, in the expressions for M25, M35, M45 and M56 make the replacement 

(1 - 2Kfe)(l - s)s ^ [(1 - 2Kb){l -p-s)-{l- 2k,){1 -p + s)]s (56) 
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The remaining four expressions read explicitly 



Mrs 
M55 

M58 



2(1 -p) -tu 



16|ks|^ + |-(1 -2ks)s + 4ks\'^ 



|l6|Kfep + 1(1 - 2Kb)s + AKbl'^ + p 

+ 16 Re {SpKbK* + s[Kb- 2(1 + p)K,bK,l + p/t*]} , 
2 f 

--Re |4(Kfe + pKs)s - [4p(/t5 + + s((l - 2/t,) + p(l - 2/tfc))] _ 
^Re |4(/t6 + p/t,)s - [4p(/tb + Ks) + s((l - 2k,) + p(l - 2Kfc))] ^^^^_ 



(57) 

1-.) 

(58) 
, (59) 



(60) 



Here, s, t and m are the kinematical variables defined in (p6|), x = x/(l — p) with x = s,t,u 
p = ml /ml and Kg = k ( j with k(s) from (p!7|). M77 and Mgs are already given in ( pSj) 



and all other contributions vanish. 

References 

[1] R. Ammar et al. (CLEO Collaboration), Phys. Rev. Lett. 71 (1993) 674. 

[2] M.S. Alamet al. (CLEO Collaboration), Phys. Rev. Lett. 74 (1995) 2885. 

[3] J. Chay, H. Georgi and B. Grinstein, Phys. Lett. B247 (1990) 399, 

I. Bigi, N. Urahsev and A. Vainshtein, Phys. Lett. B293 (1992) 430, (E) B297 (1993) 477. 
A. Falk, M. Luke and M. Savage, Phys. Rev. D49 (1994) 3367. 

[4] M. Neubert, Phys. Rev. D49 (1994) 4623. 

M. Shifman, N.G. Urahsev and A.I. Vainshtein, Int. J. Mod. Phys. A9 (1994) 2467. 

[5] R.D. Dikeman, M. Shifman and N.G. Urahsev, preprint TPI-MINN-95/9-T, UMN-TH- 
1339-95, UND-HEP-95-BIG05 ||hep-ph/9505397|l . 



[6] A. Ah and C. Greub, Phys. Lett. B361 (1995) 146. 
[7] A. Kapustin and Z. Ligeti, Phys. Lett. B355 (1995) 318. 
[8] A. Kapustin, Z. Ligeti and H.D. Politzer, Phys. Lett. B357 (1995) 653. 
[9] A. Ah and C. Greub, Z. Phys. C49 (1991) 431. 
[10] B. Grinstein, R. Springer and M.B. Wise, Nucl. Phys. B339 (1990) 269. 



21 



[11 

[12, 



[13 
[14 

[is: 

[16 

[ir; 

[18, 



A.J. Buras, M. Misiak, M. Miinz and S. Pokorski, Nucl. Phys. B424 (1994) 374. 

A. Ali, C. Greub and T. Mannel, Rare B decays in the Standard Model, in: Proc. ECFA 
Workshop on a European B Meson Factory, eds. R. Aleksan and A. Ali, DESY 90-053 
(1993). 

K. Adel and Y.-P. Yao, Phys. Rev. D49 (1994) 4945. 

M. Ciuchini, E. Franco, G. Martinelli, L. Reina and L. Silvestrini, Phys. Lett. B334 (1994) 
137. 

R. Gastmans and R. Meuldermans, Nucl. Phys. B63 (1973) 277, 
W.J. Marciano and A. Sirlin, Nucl. Phys. B88 (1975) 86. 

For a thorough discussion, see T. Muta, Foundations of Quantum Chromodynamics, World 
Scientific Publishing (1987), and references therein. 

W.A. Bardeen, A.J. Buras, D.W. Duke and T. Muta, Phys. Rev. D18 (1978) 3998. 

G. Buchalla, A.J. Buras and M.E. Lautenbacher, preprint MPI-Ph/95-104, TUM-T31- 
100/95 [|hep-ph/9512380|l , to appear in Rev. Mod. Phys. 



[19] C. Greub, T. Hurth and D. Wyler, preprint SLAC-PUB-96-7114, ZU-TH 7/1996 gep 
I ph/9603404 



22 



